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Abstract 

In an accelerating universe in General Relativity there is a maximum radius 
above which a shell of test particles cannot collapse, but is dispersed by the 
cosmic expansion. This radius could be used in conjunction with observations 
of large structures to constrain the equation of state of the universe. We 
extend the concept of turnaround radius to modihed theories of gravity for 
which the gravitational slip is non-vanishing. 
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1. Introduction 

Since 1998 cosmologists and theoretical physicists have been trying to 
explain the present acceleration of the universe discovered with type la su¬ 
pernovae. Apart from the problematic cosmol^ical constant, explanations 
based on a dark energy introduced ad hoc (see |l[ for a review) are not satis¬ 
factory and many researchers have turned to conteimlating the possibility of 
modifying gravity at large scales (j^, see BHBBB for reviews). Although 
modifying gravity is a viable possibility, too many dark energy and modihed 
gravity models ht the observational data and it is important to use any test 
of gravity which may become available, at all scales, to obtain hints on the 
correct explanation of the cosmic acceleration and, possibly, on the correct 
theory of gravity. 

One possibility which has been pointed out recently is testing theoret¬ 
ical predictions of the turnaround radius with astronomical observations 
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24 |25[ |. In an 


0, B 0 [U, H 0 Q 0 0 0 0, Eg 

accelerating Friedmann-Lemaitre-Robertson-Walker (FLRW) universe, there 
is a maximum (areal) radius beyond which a spherical shell of dust cannot 
collapse but expands forever, driven by the cosmic accelerated expansion. We 
formulate our hnal results in terms of the areal radius {R) because solutions 
of modihed gravity theories are reported in the literature using various ra¬ 
dial coordinates. However, effects peculiar to a mrticular coordinate system 
would not be meaningful in relativistic gravitypj while effects characterized 
in a geometric, coordinate-independent, way are physically meaningful. The 
areal radius separating two points in space is a physical distance identihed 
in a completely geometric way by the area of 2-spheres of symmetry in a 
spherically symmetric spacetime. In a spatially FLRW universe with line 

element = —dR + a?‘{t) {dr'^ + ^ the physical (areal) radius is 

R(f, r) = a{t)r and it expands with the scale factor a(f), while the comoving 
radius r is simply a label attached to elements of the cosmic fluid. 

The hrst comparisons of the prediction for the turnaround radius in the 
ACDM model of General Relativity with objects in the sky have been carried 
out [23|, 124, 124, 122 I and, although the precision is still poor, the method holds 


promise. In General Relativity, the concept of turnaround radius can be made 


more rigorous by using the Hawking-Hayward quasi-local mass [2a|- Given 


the motivation for modihed gravity in cosmology, here we propose to extend 
the scope of studies of the turnaround radius to alternative theories of gravity. 
We do not commit to any specihc modihed gravity theory at this stage, but 
adopt a post-Friedmannian approach in which a post-FLRW metric hts 
many theories, to lowest order in metric perturbations from an exact FLRW 
background. Gontrary to General Relativity, in which two scalar potentials 
in the perturbed FLRW metric coincide (apart from the sign), in modihed 
gravity there are two distinct potentials which are not trivially related. We 
derive the turnaround radius in this scheme and hnd that the physical (areal) 
turnaround radius depends on both potentials. We point out that, in order 
for studies of the turnaround radius to be meaningful, it is not sufficient to 
pick a theory of modihed gravity but ehorts must be made to establish which 
solutions of this theory are generic in some appropriate sense. 

The metric signature employed in this paper is —h -1--|- and we use units 
in which Newton’s constant G and the speed of light c assume the value 


’^Unless that coordinate system is associated with a preferred family of observers. 
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unity. 


2. Turnaround radius in the parametrized post-Priedmannian ap¬ 
proach 


The parametrized post-Friedmannian approach [27[ describes perturba¬ 
tions of a FLRW universe in theories of gravity alternative to General Rela¬ 
tivity. The line element ([1]) below is a rather general parametrization of the 


metric describing perturbed FLRW universes in modihed gravity j28|, l29|, [30 
(it holds, for example, in f{R) CTav ity 
conformal Newtonian gauge is |33|, [SJ, 


31 . 321). The spacetime metric in the 


ds^ = a?{'q) [— (1 -t- 2ip) drf + (1 — 20) {dr"^ + r^(iG( 2 ))] 


( 1 ) 


where + sin^ 6 dip"^ is the metric on the unit 2-sphere, rj is the con¬ 

formal time (related to the comoving time t hy dt = adrj), a{ri) is the scale 
factor of the spatially flat FLRW background, and 0 and 0 are two post- 
Friedmannian scalar potentials. While in General Relativity it is 0 = —0, 
in many modihed theories of gravity these two potentials do not coincide 
in absolute value and the gravitational slip ^ = (0 —0)/0 is used in ex- 

g iments aiming at detecting deviations from the standard AGDM scenario 
40, 35, [H, 37, ^|. The ansatz ([1]) does not include vector and ten¬ 

sor metric perturbations, which is justihed at lower order for non-relativistic 
velocities, and is common practice in the literature on cosmological pertur¬ 
bation^ {e.g., jill, ^1). 


25, 22 


^Lowing the literature on the turnaround radius in cosmology j23|, l2J, 
we assume that the FLRW perturbation is spherically symmetric. 


i.e., 0 = = 0(r). The numerical importance of deviations from 


spherical symmetry was discussed in |^, |23|, |2^. The simplest dehnition 


of turnaround radius consists of considering spherical shells of test particles 
respecting the spacetime symmetry (that is, expanding or contracting but 


^An effect due to a metric of the form m would signal modified gravity and, neglecting 
vector and tensor degrees of freedom in the metric (which is legitimate at this order of 
approximation), the potentials 0 and 0 are all that remains in the line element. However a 
better characterization of modified gravity than the metric is, at least in principle, needed 
and a rigorous derivation of the line element m in various classes of modified gravity 
theories is still missing. Nevertheless, the turnaround radius will probably be unaffected 
by these theoretical improvements. 
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not shearing nor rotating) with areal radius -R, and imposing that they have 
zero radial acceleration, R = 0, where an overdot denotes differentiation with 
respect to the comoving time t of the FLRW background. We will adopt here 
this common criterion which, in General Relativity, can be justified rigorously 

a by making use of the Hawking-Hayward quasi-local mass construct [^, 
In general, this quasi-local energy construct is not defined in modified 
gravity and here we will stick to the simple definition R = 0 for shells of test 
particles (dust) in radial motion. 

Massive test particles follow timelike geodesics with 4-tangents satis¬ 
fying UcU^ = — 1 and the geodesic equation 


du°' 

dr 


+ = 0 , 


( 2 ) 


which we choose to be affinely parametrized by the proper time r, and where 


1 


Tfoc 2^ i.9db,c ff 9dc,b 9bc,d) 


(3) 


are the coefficients of the metric connection. We will perform first order 
calculations in the metric perturbations 0 and 0, but the density contrast 
related to the spatial Laplacian of these metric perturbations is allowed to 
be large MjSlI, 421. The normalization 


= 9oo{u^Y + 9ii{u^Y = -1 


(4) 


for massive test particles with purely radial motion = 0) yields 

(«“)^ = 4 (1 - 2^.) + (1 - 2.#. - 2i,) (6) 

to first order. Eq. (|2]) then gives 

— + = o (6) 

for /i = 0,1. Eq. ([3]) provides the only non-vanishing Christoffel symbols to 
hrst order 
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where a prime denotes differentiation with respect to the radial coordinate r 
and = da/drj. Therefore, it is 


+ 2'iIj'u^u^ + — (1 — 20 — 2-0) = 0 ) (9) 

dr a a 


^ ^ ^ Q _ 

dr a 


( 10 ) 


The areal radius of the spherical spacetime (a geometric and coordinate- 
independent quantity) is 


R{t,r) = ar\/l — 20 ~ ar (1 — 0) 


to first order. Since 


and 


r = 


r = 


dR 

dt 

d^R 

dR 

dr 

dt 

d 


= {dr + ar) (1 — 0) , 

= (ar -|- 2dr + ar) (1 — 0) 

dr dr] 1 dr dr v} 

dr] dt a dr dr] au^ ’ 


u 


dt V au^ 


a 1 d d 

a? a?u^ dr 1 


we have 


d^R 

dR 


du^ 1 d f 

0 0 3“ \ ~o 

au^ aw’ dr \ w’ 


(1-0) 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


The criterion d?R/dR = 0 locating the turnaround radius j23l.l24l25| becomes 


TT “ 1 

ar + H — + 


1 dv} 


u 


dvP 


au^ dr (m^)^ dr 


= 0 , 


(17) 


where H = a/a is the Hubble parameter in comoving time. To zero order we 
have R = ar and R = dr+ ar = HR (the Hubble law) for timelike geodesics, 
so that their 4-tangents have components 


+ 5u^ = I -, 0, 0, 0 ) -|- 5u^ 


(18) 
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in coordinates {T],r,6,ip), where the perturbations 6u^ are of hrst order. To 
zero order it is = drj/dr = dr]/dt = 1/a and = 0. Eq. flT7|) now 
yields, to hrst order, 

dr + 2d5u^ + a ^ = g . (19) 

dr 

Eq. fflSj) gives = a~^ and = 6u^ which, substituted into eq. ffTOjl . yields 


d{6u^) ilj' 2a 


dr 




( 20 ) 


Inserting this equation into eq. 0191) and using the fact that = ad hnally 
gives 

( 21 ) 


ar 


lb' 

— = 0 


to hrst order. Eq. fl?I|) locates the comoving turnaround radius Vc once the 
post-Friedmannian potential ip is given. In general, this is a trascendental 
equatiorj^ and reduces to an algebraic one only for simple forms of ip. In an 
accelerating universe (a > 0), and assuming a positive decreasing function 
ip'{r) (as is the case, for example, for a Newtonian point mass ip{r) = —m/r 
and as it is reasonable to expect in general for a single spherical perturbation), 
there exists a unique intersection between the straight line of positive slope 
y = dr and the decreasing function y = ip'{r), which dehnes uniquely the 
turnaround radius re- In a decelerating universe (a < 0), by contrast, there is 
no intersection between these two curves and there is no turnaround radius. 

Eq. fl2T|) is the main result of this paper. The job is not complete, though, 
because one would like to know the areal, not the comoving, turnaround 
radius. The areal turnaround radius is 


Rc = a{t)rc [1 - (p{rc)] , ( 22 ) 

where Tc is the unique root of eq. fl^T]) . The comoving turnaround radius 
depends only on the post-Friedmannian potential ip but the areal turnaround 


^ Often in Brans-Dicke and f{R) gravity one finds TZi iM 13| potentials ip oc {y/rY 
where /r is a constant mass coefficient and the exponent p is not integer, which makes 
eq. (ICTl non-algebraic. 


6 









radius, which is the physical radius, depends also on the second potential 0. 
By expanding this equation as 


R 

a 




(23) 


to hrst order, eq. fl^ turns into the equation satished by the areal turnaround 
radius Rc 

dRc + dRc(j)c ~ ~ 0 ( 24 ) 

where 0c = 4>{Rc/ci) and 0' = 'ip'{Rc/a), which makes it clear that Rc 
depends on both potentials 0 and 0. To make explicit the difference between 
modified gravity and General Relativity one can use the gravitational slip 
which relates the potentials 0 and 0 through 0 = 0(1 — ^). Eq. fl2T|) for the 
turnaround radius Vc in the geometry ([T]) can be written as 


arc 


a a 


0 . 


Correspondingly, eq. fl2Tll becomes 


(25) 


iiR,{l + 4>c)-4>'c{i-(c)+<l‘A = 0- (26) 


3. Examples 

Here we consider a few examples. Unfortunately, only a handful of ana¬ 
lytical solutions of modified gravity theories are known which describe spheri¬ 
cally symmetric inhomogeneities embedded in FLRW universes, and for these 
spacetimes the cosmological “background” is invariably decelerated instead 
of accelerated (see the recent review j^), so they are not useful here. 

3.1. ACDM model 

As a first example, consider the ACDM model of General Relativity and 
the simple spherical perturbation potentials 0(r) = —4>{r) = —m/r, with m 
a mass constant. This expression makes sense for sufficiently large values of r 
in order to keep the metric perturbations small. Using the FLRW acceleration 
equation 

- = -0 (3F + p) , (27) 
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where P and p are the pressure and energy density of the dark effective fluid 
dominating the dynamics of the accelerated universe, respectively, eq. fET]) 
provides the comoving turnaround radius 


Tr = 


3m 


47r Ip + 3P\ 


1/3 


(28) 


The corresponding areal radius is, to dominant order, 

3ma \ 

^Att Ip + 3P\) 


Rr = avn = 


(29) 


where the comoving scale ma, rather than the constant m, is the physical 
mass to consider (see the corresponding mass in Refs. 23|, |2J, |2^ and the 
discussion in ji^). This equation reproduces previous expressions of the 
turnaround radius found in Refs, j^, SH, EH, El, 13, 14, 15, 16, 17, 18, 19 


20l. l2ll. I 22 I. I 23 I. I 24 , l25|. In an expanding universe, the quantity ma is identified 


with the total mass enclosed in a sphere of comoving radius r and areal radius 
R = ar. This mass comprises the contribution of the “local” perturbation 
plus that of the cosmic fluid enclosed in the sphere IH, Hi . In General 
Relativity, this mass ma coincides with the Hawking-Hayward quasi-local 
energy |H, l46| and, because of spherical symmetry, also with the Misner- 
Sharp-Hernandez mass H|. In modihed gravity, however, these quasi-local 


energy constructs are not defined. 

3.2. Generalized McVittie attractor 

Generalized McVittie solutions of the Einstein equations were introduced 
The old McVittie class of solutions of General Relativity 


m 


51 


53 de¬ 


scribes a central object (under certain conditions, a black hole) embedded in 
a FLRW universe. The McVittie class of spacetimes has been the subject of 
much recent attention [H, 55, [H, IH, 59, 0, 611 and it has been shown 
to solve the equations of cuscuton theory, a special Hofava-Lifschitz theory 


62j . Accretion of the surrounding cosmic fluid onto the central McVittie ob¬ 
ject is artificially forbidden: in McVittie imposed that the time-radius 
component Gtr of the Einstein tensor vanishes, corresponding to zero radial 
flux of mass-energy onto the central inhomogeneity. The rationale for this 
assumption was simplifying the search for a solution of the Einstein equations 
describing a central object (in McVittie’s intentions, a point mass) embedded 
in a cosmological space. McVittie realized that it was much simpler to forbid 


























































accretion than having to model it. Nonetheless, in General Relativity, one 
can allow a radial flow of energy, thus obtaining the so-called “generalized 
McVittie solutions” j^. This radial energy flow is spacelike and rather ar¬ 
tificial in this context. All simple models of heat conduction or of viscosity 
of imperfect fluids in relativity suffer from the inconsistency that heat con¬ 
duction is instantaneous and obviously contradicts relativity, but imperfect 
fluids are used routinely as simple models because a proper relativistic ther¬ 
modynamical treatment increases the complication to the point of rendering 
simple models impossible. Within the context of Horndeski gravity, the gen¬ 
eralized McVittie solutions become more realistic scalar held solutions |63|, 
in which the scalar held is suceptible of an effective fluid description. Within 
the class of generalized McVittie solutions there is a late-time attractor 64 
with geometry described by 


ds^ = 


2mr 


dt^ -I- a^{t) 


dr"^ 


_ 2mQ 


+ r^dn^ 


( 2 ) 


(30) 


where mo is a constant. In the approximation mo/r 1, this line element 
reduces to 


ds^ = 


2mr 


dt^ -|- a^{t) 


l + ^\dr‘^ + r^dn^ 


( 2 ) 


(31) 


with -0 = —m/r = —cj) and we reduce formally to the previous example. This 
relation is not surprising since the metric fl30|) solves both the Einstein and 
the Horndeski held equations and the ansatz ([1]) represents a solution of the 
held equations of a theory, not the theory itself. Therefore, in this case, the 
turnaround radius does not discriminate between the ACDM model based 
on General Relativity and Horndeski gravity. 


4. Discussion and conclusions 


In generalized gravity, when the line element assumes the post-Friedman- 
nian form ([T]) widely used in the literature and with a presently accelerating 
FLRW background, we have found a formula (eq. (El])) for the turnaround 
radius Re- The upper bound to the radius of stable virialized structures 
that can be obtained by observing the largest bound objects in the universe 
(those that are on the verge of breaking apart) would constitute an observable 
to look at in order to see a signature. Indeed the potentials (j) and "0 are 


commonly regarded as observables in modified gravity [6a, 1^ 
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Several modified gravity theories admit multiple representations, the most 
well known situation being that of scalar-tensor gravity which can be ana¬ 
lyzed in two conformal frames, the so-called Jordan frame and Einstein frame 
related by a conformal transformation of the metric plus a non-linear redef- 
inition of the Brans-Dicke-like scalar held present in the theory [^, ^ 


69 


These two frames are in principle physically equivalent provided that a rescal¬ 
ing of the units of mass, length, and time (and of all the derived units) is 
performed 


see also ITO 


7ll | and the references therein). However, it may 
be rather cumbersome to relate physical quantities in the two frames (see, 
e.g., the recent discussion [72|), which renders the equivalence less useful for 
practical purposes. 

Before eq. can be compared with data on observed structures which 


are presumably on the verge of virialization (as done in Refs. 23|, |2J, |25|). 
some more work needs to be done. Eq. fl2T]) may test how accurately eq. o 
describes the spacetime geometry to hrst order in the metric perturbations, 
and it is unlikely that observations will require higher precision than hrst 
order. However, the same spacetime metric may be a solution of several 
theories of gravity, including General Relativity, at once. Thus, the the¬ 
oretical problem arises of which solutions of a given theory of gravity are 
signihcant, in the sense of being generic and of being different from the 
standard one of Einstein’s theory. There is some degeneracy, in the sense 
that different theories of gravity, including General Relativity, can admit the 
same post-Friedmannian metric as a solution to hrst order in the metric per¬ 
turbations (the same way that most theories of gravity admit the FLRW 
or the Schwarzschild metric as an exact solution). If, in addition to the 
FLRW background, the perturbed solution is the same then observational 
constraints on the turnaround radius will not detect deviations from Gen¬ 
eral Relativity even if gravity was described by an alternative theory. Work 
needs to be done to understand which perturbed cosmological solutions are 
in some sense “generic” and, therefore, expected in various theories of grav¬ 
ity and whether the turnaround radius computed for such theories dihers 
signihcantly from that of Einstein’s theory. If modihed gravity does indeed 
explain the cosmic acceleration without dark energy, as is hypotesized, then 
it deviates from Einstein gravity at large scales and the turnaround radius 
holds promise for an independent test of gravity at the interface between 
astrophysics and cosmology. 

A completely different, statistical mechanics approach to the turnaround 
radius is contained in Refs. [73|, in which galaxy clusters are studied by 
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including effective “dark energy particles” forming a perfect fluid and it is 
found that there is a “reentrant radius” corresponding to our turnaround 
radius. Beyond this radius thermal equilibrium of the gas is not possible 
and quintessence (with P > —p) increases the gravitational instability, while 
phantom energy (with P < —p) helps stabilizing against gravitational col¬ 
lapse. A comparison betwen our macroscopic approach and the statistical 
mechanics approach of j^, 74] is not straightforward because these authors 
assume a barotropic equation of state of state for dark energy which is linear 
and constant. While it is possible to write the held equations of modihed 
gravity in the form of effective Einstein equations by collecting geometric 
terms into the right hand side as an effective huid, in general the latter is not 
a perfect huid and the ehective equation of state is not linear nor constant. 
Imposing a constant ehective equation of state would results in unphysical 
restrictions on the theory of gravity motivated only by mathematics. Never¬ 
theless, it will be interesting to generalize the approach of 73|, |7^ to modihed 
gravity in the future. 
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